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ABSTRACT 

We studied the Active Galactic Nuclei (AGN) radio emission from a compilation of hard X-ray 
selected samples, all observed in the 1.4 GHz band. A total of more than 1600 AGN with 2-10 keV 
de-absorbed luminosities higher than 10**^ erg s~^ cm~^ were used. 

For a sub-sample of about 50 z<Q.l AGN it was possible to reach a ~80% fraction of radio detec- 
tions and therefore, for the first time, it was possible to almost completely measure the probability 
distribution function of the ratio between the radio and the X-ray luminosity Y{,x= log(Li.4/Lx), 
where Li.4/Lx= vL^ilA GHz)/Ljf (2-10 keV). The probability distribution function of Kx was func- 
tionally fitted as dependent on the X-ray luminosity and redshift, P{Rx\Lx, z). It roughly spans over 
6 decades (— 7<Rx< — 1), and does not show any sign of bi-modality. It resulted that the probability 
of finding large values of the Kx ratio increases with decreasing X-ray luminosities and (possibly) 
with increasing redshift. No statistical significant difference was found between the radio properties 
of the X-ray absorbed (Nh>10^^ cm~^) and unabsorbed AGN. 

The measure of the probability distribution function of Rx allowed us to compute the kinetic 
luminosity function and the kinetic energy density which, at variance with what assumed in many 
galaxy evolution models, is observed to decrease of about a factor of five at redshift below 0.5. About 
half of the kinetic energy density results to be produced by the more radio quiet (Rx<-4) AGN. 

In agreement with previous estimates, the AGN efficiency Ckin in converting the accreted mass 
energy into kinetic power (L;<-=efci„mc^) is, on average, ekin—^ x 10~^. The data suggest a possible 
increase of ekin at low redshifts. 

Subject headings: galaxies: active - galaxies: evolution - galaxies: luminosity function, mass function 
- quasars: general - radio continuum: galaxies - X-rays: galaxies 



1. INTRODUCTION 

One of the main questions in galaxy formation evo- 
lution studies is the role of AGN feedback. Accord- 
ing to popular AGN/galaxy co-evolutionary scenarios, 
I once central super massive black holes (SMBHs) reaches 
. masses >10'^ - 10^ Mq, the AGN can heat efficiently the 
galaxy interstellar matter (ISM) through winds, shocks, 
and high-energy radiation (Silk & Rees 1998, Fabian 
1999), inhibiting further accretion and star-formation 
' and making the galaxy colors redder (see e.g. Cattaneo 
et al. 2009 for a review). 

Unfortunately there are still relatively few direct ob- 
servations of AGN feedback, and its inclusion in galaxy 
evolution models is often performed using adjustable 
parameters to obtain the observed galaxy colors. In- 
deed, the results of the hydrodynamic N-body simula- 
tions (see, e.g. Di Matteo, Springel & Hernquist 2005; 
Springel 2005; Hopkins et al. 2005, 2006) and of the 
semi-analytical models of galaxy formation and evolu- 
tion (SAMs, Monaco, Salucci & Danese 2000; Kauffmann 
& Haenhelt 2000; Volonteri, Haardt & Madau 2003; 
Granato et al. 2004; Menci et al. 2006, 2008; Croton et 
al. 2006; Bower et al. 2006; MarulH et al. 2008) depend 
on the AGN triggering mechanism and on the AGN feed- 
back description. Is the feedback at work mainly during 
luminous AGN phases (the so called "quasar" or "radia- 
tive" mode, see e.g. Menci et al. 2008)? Or rather AGN 



feedback proceeds continuously during the Cosmic time 
at a low rate (the so called "radio" or "kinetic" mode, 
see e.g. Croton et al. 2006, MaruUi et al. 2008)? 

In the first case the AGN feedback is associated to the 
the main radiative, and then X-ray, activity of the AGN 
which is thought to occur during episodes of accretion of 
cold gas coming from the galaxy ISM following galaxy 
encounters (e.g. Barnes & Hernquist 1992, Cavaliere 
& Vittorini 2000, Menci et al. 2008). In this scenario 
the AGN feedback is usually associated to the radiation 
and its total efficiency must be proportional to the AGN 
fraction (the AGN luminosity function vs. the galaxy lu- 
minosity function) and to the efficiency in releasing the 
AGN power in the galaxy ISM. 

In the second case the radio feedback is assumed to 
be related to a low, uninterrupted and constant matter 
accretion rate onto the central SMBH (~10~^ Mq/jt), 
coming from a quiescent inflow of gas cooling from the 
halo's hot atmosphere (e.g. Monaco, Salucci & Danese 
2000; Croton et al. 2006; Bower et al. 2006). This 
accretion rate is too small to contribute significantly to 
the bolometric output of the AGN. At variance, the me- 
chanical effects of accelerated particles (jets), which are 
observed to be responsible for the large cavities on the 
intra-cluster medium revealed in the X-rays (e.g. McNa- 
mara et al. 2000), are believed to significantly perturb 
the ISM into radio galaxies (Saxton et al. 2005, Suther- 



2 



LA FRANCA, MELINI & FIORE 



land & Bicknell 2007, Tortora et al. 2009; Krause & 
Gaibler 2009 and references therein). In this scenario, 
radio mode feedback total efficiency is thought to be 
proportional to the total accreted mass and then to the 
SMBH mass function and to the way the energy is chan- 
nelled by the SMBH into its host galaxy and it is released 
in the ISM and intergalactic medium (IGM). 

In the Croton et al. (2006) and Marulh et al. (2008) 
models, the radio-mode feedback is more effective in sup- 
pressing the cooling flows in the massive galaxies at late 
times (low redshifts). In the Cattaneo et al. (2006) 
model the cooling and star-formation are efficiently sup- 
pressed by the AGN radio feedback for haloes above a 
critical mass of ~10^^ Mq below z~3. Bower et al. 
(2006) assume that the AGN energy injection is deter- 
mined by a self-regulating feedback loop that starts when 
the luminosity exceeds some fraction of the Eddington 
luminosity. Recently, Shabala & Alexander (2009) have 
presented a galaxy formation and evolution model where 
the radio feedback occurs when the AGN accretion rate 
falls below a certain value and enters the advection dom- 
inated accretion flow (ADAF) regime. In their model 
the radio (and then kinetic) power is assumed to scale 
linearly with the accretion rate. 

To support the above studies several authors have com- 
pared the mechanical (kinetic) luminosity function (LF) 
of radio sources to the bolometric AGN LFs (e.g. Best 
et al. 2006; Merloni & Heinz 2008; Shankar et al. 2008; 
Kording, Jester & Fender 2008; Cattaneo & Best 2009; 
Smolcic et al. 2009). All these works are based on the 
convolution of some empirical relation between the AGN 
kinetic power and its radio luminosity (see e.g. Willott et 
al. 1999; Best et al. 2006; Merloni & Heinz 2007; Birzan 
et al. 2008) with the AGN radio LF. As a consequence, 
these works deal with the AGN as a population and do 
not allow the inclusion in the models of the kinetic power 
of each single source during its evolution. 

In the Croton et al. (2006) and Bower et al. (2006) 
SAMs it is not implemented any direct relation between 
the radio activity (feedback) and the X-ray luminosity 
(main episodes of mass accretion). However, a strong 
correlation between the X-ray (Lx) and radio luminosity 
of AGN is observed (e.g. Brinkmann et al. 2000), and 
Lx is a good proxy of the SMBH accretion rate m, via 
the knowledge of the X-ray bolometric correction Kx 
and the efficiency e of conversion of mass accretion into 
radiation 

Lboi erhc^ 
Kx (1 - t)Kx 

where L^oi is the bolometric luminosity and it is assumed 
e ~ 0.1 (see e.g. Marconi et al. 2004; Vasudevan & 
Fabian 2009). 

The aim of this paper is therefore to estimate the AGN 
kinetic power linking the AGN radio emission to the ac- 
cretion rate related to the AGN activity (the luminous 
phase). Our approach is then different from most of 
the above described models, and it can be considered, 
in some sort, as a quasar mode feedback (i.e. related 
to the luminous-accreting phases), but associated to the 
emitted radio-kinetic power. This can provide a robust 
quantitative root to a different kind of radio feedback, 
thus guiding its self-consistent inclusion in SAMs. In 
this framework, a very useful ingredient is the measure of 
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Fig. 1. — Distribution of the off-sets between the cross-correlated 
radio and X-ray sources. 

the probability distribution function P(Rx) of the ratio, 
Rx, between the AGN radio and hard X-ray luminosity 
[Rx=i^L,,(1.4 GHz)/Lx(2-10 keV)]. To this purpose we 
used a data-set of more than 1600 X-ray selected AGN, 
observed in the radio band at 1.4 GHz, to measure P(Rx) 
as a function of both luminosity and redshift: P(Rx| Lx, 
z). 

To estimate the AGN kinetic LF and its evolution, we 
first computed the radio LF by convolving P(Rx| Lx, z), 
with the AGN 2-10 keV LF. We then convolved this radio 
LF with some of the relations (available from the liter- 
ature) between the AGN radio and kinetic luminosities. 
As sanity tests, we first compared our results with previ- 
ous studies on the relationship between the AGN X-ray 
and radio luminosities, and then we checked if previous 
measures of the radio LF and counts were correctly re- 
produced. 

We adopted a fiat cosmology with Hq = 70 km 
Mpc""'^, r2j\/=0.30 and f2A=0.70. Unless otherwise 
stated, uncertainties are quoted at the 68% (la) confi- 
dence level. 

2. THE DATA 

As our objective was to use the Rx distribution in 
order to estimate the kinetic (radio) luminosity of the 
X-ray selected AGN, we had to measure a radio emission 
which was as much as possible causally linked (contem- 
porary) with the observed X-ray activity (accretion). We 
then decided to measure the radio fiuxes in a region as 
much close as possible to the AGN, therefore minimizing 
the contribution of objects like the radio lobes in FRII 
sources (Fanaroff & Riley 1974). For this reason we built 
up a large data-set of X-ray selected AGN (where red- 
shift and Nff column densities estimates were available) 
observed at 1.4 GHz with a '^1" typical spatial resolu- 
tion (in our cosmology 1" corresponds, at maximum, to 
about 8 Kpc at z~2). The cross correlation of the X- 
ray and radio catalogues was carried out inside a region 
with 5" of radius (almost less than or equal to the size of 
the central part of a galaxy like ours), following a max- 
imum likelihood algorithm as described by Sutherland 
and Saunders (1992) and Ciliegi et al. (2003). In Fig- 
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ure 1 we show the off-sets between the X-ray and radio 
positions of the whole sample. These resulted to have 
a root mean square (rms) of 1.4". Therefore we expect 
to have (properly) preferentially included compact FRI 
radio sources in our cross-correlation and have excluded 
most of the contribution of the (especially bright) ex- 
tended FRII radio lobes from our analysis (see sect 3.4.2 
for a discussion on the contribution of the excluded FRII 
sources to the radio counts) . 

The total AGN sample was built up from a compilation 
of complete (i.e. with almost all redshift and N^^ mea- 
sures available) hard (mostly 2-10 keV) X-ray selected 
samples, with unabsorbed 2-10 keV luminosities higher 
than 10^^ erg/s^, as explained below. 

2.1. The bright sample. The SWIFT, INTEGRAL and 
HEAO-1 missions 

In order to build up a large unbiased bright AGN (low 
redshift) sample, we joined the AGN samples recently 
generated from the SWIFT and INTEGRAL missions 
with the sample published by Grossan (1992) using the 
HEAO-1 data. In the case of sources in common, prior- 
ity was given first to the SWIFT data and then to the 
INTEGRAL ones and lastly to the sample of Grossan 
(1992). 

The radio information was obtained via cross- 
correlation with the 1.4 GHz radio data taken from the 
FIRST VLA survey (Becker et al. 1995). The FIRST 
images have 1.8" large pixels, a typical rms sensitiv- 
ity of 0.15 mJy and a resolution of 5". In the case of 
no radio detection a 5(7 upper limit of 0.75 mJy was 
adopted. The large area NVSS (Condon et al. 1998) and 
SUMSS/ATCA (Mauch et al. 2003) radio surveys have 
positional uncertainties significantly larger than FIRST 
and therefore were not used in our analysis. 

• SWIFT. The 5 H^/FT sample we used is composed 

by 121 sources with high galactic latitude (|6|>15°) 
and detected with 14-195 keV fluxes brighter than 
IQ-" erg s-i cm^2 (Tueller et al. 2008). All but 
one of the 121 sources have a redshift and an opti- 
cal spectroscopic classification available (Tueller et 
al. 2008). 44 out of the 121 sources have been ob- 
served by FIRST. We limited our analysis to the 40 
sources with 2-10 keV unabsorbed X-ray luminosity 
higher than 10^^ erg/s: there are 18 broad optical 
emission line type 1 and 1.5 AGN (AGNl/1.5), 13 
narrow emission line type 2 AGN (AGN2), 2 galax- 
ies and 7 BL Lac in this sub-sample. In this and 
in all the other samples, all the BL Lac were ex- 
cluded from our analysis as their radio fluxes are 
strongly amplified by the boosting of the relativis- 
tic radio jets. Therefore the AGN sample we used 
is composed by 33 sources in total (28 detected by 
FIRST). N// column densities measures were pro- 
vided by Tueller et al. (2008). 

• INTEGRAL. We tried to complement the SWIFT 
data with the catalogue of 46 sources detected by 

^ Throughout this work we assumed that all the X-ray sources 
having 2-10 keV unabsorbed luminosities higher than erg/s 
are AGN (see e.g. Ranalli et al. 2003 for a study of the typical 
X-ray luminosities of star forming galaxies). 



INTEGRAL at a 5a" significancy level by Beck- 
mann et al. (2006). However, after removing two 
sources without N^ measurements (from Sazonov 
et al. 2007), and 21 sources already included in the 
SWIFT sample, we ended with 18 sources which 
were not covered by the FIRST radio observations. 
Therefore, no source from the INTEGRAL cata- 
logue from Beckmann et al. (2006) was included in 
our analysis. 

• GROSSAN. The bright sample was eventually com- 
plemented with the sample of the HEAO-l sources 
described by Grossan (1992) as revised by Bru- 
sadin (2003). Brusadin et al. (2003) investigated, 
from the total sample of Grossan (1992), those 74 
sources with 2-10 keV fluxes brighter than 2 x 10~^^ 
erg s~^ cm~^ . For 52 of these 74 sources the origi- 
nal optical counterparts were observed in the hard 
X-rays by at least one among the ASCA, Beppo- 
SAX and XMM-Newton satellites. All the coun- 
terparts resulted to be real hard X-ray sources and 
a new estimate of the Njj column densities were 
derived. We used in this work a sub-sample of 66 
objects, from the sample of Brusadin (2003), for 
which reliable measure of the Njj column densities 
were available. 

Forty-one out of these 66 sources are not included 
in the SWIFT sample, and ten (10 AGNl/1.5) of 
them were observed by FIRST and have 2-10 keV 
X-ray unabsorbed luminosities higher than 10^^ 
erg/ s. Nine out of these last ten sources were also 
detected in the radio band by FIRST. 

In summary, the bright sample contains 43 X-ray 
sources (28 AGNl, 13 AGN2, 2 galaxies; once the BL 
Lac were excluded) observed in the 1.4 GHz radio band 
by FIRST, with 37 detections. 

2.2. AMSS 

We selected those 87 sources of the ASCA Medium 
Sensitivity Survey (AMSS; as in Akiyama et al. 2003) 
having a S/N>5.5 significant detection in the X-rays, 
|6|>30° and S>-2Q°; 76 of these sources were identified 
as AGN: seven clusters, 3 BL Lac and one star were 
excluded. Forty-three of these AGN were observed by 
FIRST (33 AGNl, 10 AGN2), and 11 were detected. 

2.3. HBSS 

We used the 67 sources of the XMM-Newton Hard 
Bright Sensitivity Survey (HBSS; Delia Ceca et al. 
2004) with 4.5-7.5 keV ifluxes brighter than 7x10'^'^ 
erg cm~^ (which corresponds to a 2-10 keV limit 
of 3.5x10"^^ erg s~^ cm~^ if a spectral index a = 0.7, 
where Fi, oc is assumed). A subsample of 62 sources 
were used in our analysis after the exclusion of 2 stars, 
one cluster and 2 sources without a spectroscopic iden- 
tification. Thirty-two sources were observed by FIRST 
(23 AGNl, 9 AGN2), while 6 were detected. 

2.4. ALSS 

The ASGA Large Sky Survey (ALSS) is a contiguous 7 
deg^ strip in the North Galactic Pole region (Ueda et al. 
1999); we selected a sample with a limiting 2-10 keV flux 



4 



LA FRANCA, MELINI & FIORE 



TABLE 1 
The Samples 



Sample 


Nl.4 


Nx 


«1.4 
Nx 


Fl.4 


Fx 




i-'l.i 
Fx 




(1) 


(2) 


(3) 


(4) 


(5) 




(6) 


Bright 


37 


43 


0.86 


750 


2x10" 


-11 


-6.3 


AMSS 


11 


43 


0.26 


750 


3x10" 


-13 


-4.5 


ALSS 


9 


30 


0.30 


750 


1x10" 


-13 


-4.0 


BBSS 


6 


32 


0.19 


750 


3.5x10" 


-13 


-4.5 


COSMOS 


141 


677 


0.21 


75 


3x10" 


-15 


-3.5 


CLANS 


69 


139 


0.50 


19 


3x10" 


-15 


-4.1 


ELAIS-Sl 


45 


421 


0.11 


150 


2x10" 


-15 


-3.0 


CDF-S 


12 


94 


0.13 


70 


2.6x10" 


-16 


-2.4 


CDF-N 


45 


162 


0.28 


45 


1.4x10" 


-16 


-2.3 


Total 


375 


1641 


0.23 











Notes. 1) Number of AGN detected at 1.4 GBz; 2) Number of 
X-ray AGN observed at 1.4 GBz; 3) Fraction of AGN detected at 
1.4 GBz; 5) 1.4 GBz radio flux limit in /ijy units; 5) 2-10 keV flux 
limit in erg s"! cm-^ units; 6) log[!^F,41.4 GBz)/F^{2-10 keV)]. 
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Fig. 2. — Lower (5it) radio flux 1.4 GBz limits as a function 
of the lower X-ray 2-10 keV flux limits of the samples. Dashed 
lines show the loci with equal i^F^(1.4 GBz)/F^(2-10 keV) ratio as 
shown by the labels. 



of 1x10^^^ erg s"^ cm"^ from Akiyama et al. (2000). 
This sample contains 30 AGN (25 AGNl, 5 AGN2), as 
well as two clusters, one star and one object without 
spectroscopic identification, which were excluded from 
our analysis. All the 30 AGN were observed by FIRST, 
while nine were detected. 

2.5. COSMOS 

In order to use the COSMOS survey, we used the 
XMM-Newton X-ray catalogue by Cappelluti et al. 
(2009) with a limiting 2-10 keV flux of -3x10-1^ 
erg s~^ cm~^ and cross-correlated with the spectroscopic 
identifications by Trump et al. (2009) and the photo- 
metric redshift estimates from Ilbert et al. (2009) and 
Salvato et al. (2009). The Nh measures were derived 
by our analysis of the hardness ratios (HR). The radio 
data were obtained from Schinncrcr et al. (2007). In or- 
der to allow a uniform radio coverage (with 1.4 GHz rms 
of 15 /iJy) the central 1 deg^ squared area with limits 
9''58'40"<a<10''2'40" and 1°42'<(5<2°42' was used. In 
this area the X-ray catalogue contains 712 sources. 677 
are extragalactic with a redshift measure available and 
2-10 keV luminosity higher than 10**^ erg/s. 389 have a 
spectroscopic redshift (186 AGNl, 52 AGN2, 71 emission 
line galaxies, ELG, 32 normal passive galaxies, GAL, 48 
no class) , while 288 have only a photometric redshift esti- 
mate available. We used all the radio detections with 1.4 
GHz fiux limits brighter than 75 fiJy (5cr) and used the 
same threshold as an upper limit for all the remaining 
sources even though lower flux detections were available 
in some cases. In total 141 out of the 677 sources, con- 
tained in our selected region of the COSMOS field, were 
detected in the 1.4 GHz radio band. 

2.6. ELAIS-Sl 

In the European Large Area ISO Survey field SI 
(ELAIS-Sl) we used the catalogue of the XMM-Newton 
sources published by Puccetti et al. (2006), which 
reaches a 2-10 keV fiux limit of 2x10^^*^ erg s~^ cm~^. 
The spectroscopic identifications and classifications pro- 
vided by Feruglio et al. (2008) and Sacchi et al. (2009) 
were used, while the 1.4 GHz radio data (with a 5ct limit 
of 150 /iJy) were taken from Middelberg et al. (2009). 




2 3 
redshift 

Fig. 3. — 2-10 keV X-ray de-absorbed luminosity as a function of 
the redshift for all the AGN observed at 1.4 GBz used in this work. 
Top. Crosses indicate the sources observed in the radio band, while 
flUed squares show those objects which have been detected in the 
radio. The dashed lines show the two z-hx regions used to measure 
the distribution of Hx at Rx<-4. Bottom. Crosses indicate the 
sources with log N^f < 22 cm~^, while filled squares show those 
AGN with log Nh > 22 cm~^. The dashed lines show the five 
z-hx regions used to measure the distribution of Rx at 

The whole sample contains 421 extragalactic sources 
with a redshift measure available and 2-10 keV lumi- 
nosity higher than 10^^ erg/s: 240 have been identified 
and classified spectroscopically (116 AGNl, 34 AGN2, 
68 ELG, 22 GAL), while 181 have a photometric red- 
shift available. Forty-five out of these 421 AGN were 
detected at 1.4 GHz above the 150 /iJy limit. 
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2.7. CLANS 

In the Chandra Lockman Area North Survey (CLANS) 

field the X-ray data were taken from Trouille et al. 
(2008), which also publish the spectroscopic and pho- 
tometric rcdshift measures. Our sample consists of the 
sources with S/N>3 detections in the hard 2-8 keV band 
and includes in the circular area with radius of 0.32 deg 
and center in: a 10''46', 5 59°01' (J2000). This area con- 
tains 139 extragalactic sources with a redshift measure 
available and 2-10 kcV luminosity higher than 10^^ erg/s: 
113 sources were spectroscopically identified (58 AGNl, 
34 AGN2, 16 ELG, 5 galaxies), while 29 sources have a 
photometric redshift available. The Nij measures were 
derived from the analysis of the HR. Radio data were 
obtained from Owen ct al. (2008). We conservatively 
modeled the spatial dependence of the (Str) radio flux 
limits, which vary from 18.5 jiZy in the central region up 
to 59.2 /iJy near the edges. Sixty-nine out of the 139 
sources were radio-detected. 

2.8. CDFS 

In the Chandra Deep Field South (CDFS) the subsam- 
ple of the GOODS-S X-ray sources from the catalogue 
of Alexander et al. (2003) was used. The whole sample 
consists of the 94 point like extragalactic sources over 
a 2-10 kcV flux limit (at the aim point) of 2.6 xlO"!*' 
erg s~^ cm^^ with a redshift measure available and a 2- 
10 keV luminosity higher than 10''^ erg/s. Spectroscopic 
identifications for 69 sources, as well as 25 photomet- 
ric redshifts, were obtained from (Brusa et al. 2009b). 
Twenty-nine sources were idcntifled as AGN (17 AGNl, 
12 AGN2), 39 as ELG and one as a galaxy. The radio 
data were taken from Miller et al. (2008) and have a 
1.4 GHz 5(7 flux limit of 70 /zJy. Twelve sources were 
radio-detected. 

2.9. CDFN 

In the Chandra Deep Field North (CDFN) we used the 
296 2-8 kcV X-ray sources detected by Chandra within 
the GOODS-N area with a 2-10 keV flux hmit (at the 
aim point) of 1.4x10""'^^ erg s~^ cm^^ by Alexander ct 
al. (2003). To convert the 2-8 keV fluxes in the 2-10 keV 
band a spectral index a=0.4 were assumed. The spec- 
troscopic idcntiflcations were taken from Trouille et al. 
(2008). The sample consists of 162 extragalactic sources 
with a redshift measure available and 2-10 keV luminos- 
ity higher than 10"^^ erg/s: 104 sources were spectro- 
scopically identified (16 AGNl, 20 AGN2, 46 ELG, 22 
galaxies), while 50 sources have a photometric redshift; 
eight other spectroscopical redshift were retrieved from 
literature. N// values were obtained from the analysis of 
the HR. The radio information were obtained from the 
new data reduction from Biggs & Ivison (2006) of the 1.4 
GHz VLA observation of Richards (2000). Forty- five out 
of the 162 extragalactic sources were identified in the 1.4 
GHz band over a (5cr) flux limit of 45 Jy. 

2.10. The whole X-ray sample 

In summary, our radio observed hard X-ray selected ex- 
tragalactic total sample contains 1641 sources with both 
redshifts (either spectroscopic or photometric) and Njy 
column densities measured, and with un-absorbed 2-10 
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Fig. 4. — Probability distribution function of Kx for the full 
sample under the assumption of a flat distribution (dashed line). 

keV luminosities higher than 10"*^ crg/s^. 1003 sources 
have N/f higher than 10^^ cm^^ (hereafter deflned "X- 
ray absorbed"). 375 (23%) sources were detected in the 
1.4 GHz radio band. See Table 1 for a summary of the 
main properties of all samples used. In Figure 2 we show 
the (5(t) radio 1.4 GHx flux limits of each survey as 
a function of their deepest 2-10 kcV X-ray flux limits. 
The distribution of the de-absorbed X-ray luminosity of 
all sources (distinguished according to radio detection 
and X-ray absorption classes) as a function of redshift is 
shown in Figure 3. 

3. THE DISTRIBUTION OF THE R=Lfl/Lx RATIO 

3.1. The method 

We searched for a functional fit of the probability dis- 
tribution function of Rjc , as a function of the X-ray lumi- 
nosity, Lx, and the redshift: P{Rx\Lx,z). The method 
is based on the comparison, through estimators, of the 
observed and expected numbers of AGNs (in the Lx-z- 
Rjf space) obtained by taking into account the obser- 
vational selection effects (i.e. the radio flux limits) of 
each sample. Once a probability distribution function 
P{Rx\Lxtz) is assumed, the number of expected AGNs 
in a given bin of the Ljsf-z-Rx space is the result of the 
sum, over the number of all AGN contained in the hx-z 
bin, of the expected number of AGNs contained in that 
bin of Rx, by taking into account the radio flux limits 
on each source. This method reproduces the observa- 
tions and consequently properly takes into account both 
the radio detections and the upper limits (see La Franca 
et al. 1994, 1997, 2002, 2005 for similar applications). 

3.2. The fit 

As a first test we assumed a constant (fiat) probability 
distribution function of Rx in the range -7<Rjjc<0 (see 
Figure 4). This distribution, although different from the 
true one, allows to see which is the shape of the aver- 
age true distribution via the analysis, in each bin, of the 
deviations of the observed numbers of AGN from the 

2 Other X-ray samples, sueh as HELLAS2XMM (Fiore et al. 
2003; Cocchia et al. 2007), the XMM/Lockman Hole (Brunner et 
al. 2008) and the XMM Medium Survey (Barcons et al. 2007), 
were not included in our work because unbiased, homogeneous and 
dedicated radio 1.4 GHz observations were not available. 
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TABLE 2 
Best fit parameters 



Model 


N 


7h 


7l 


Rc 


OIL 




xVd.o.f. 




1 - box 














670.17/12 





2 - no dep 


1.0620 


0.476 


1.93 


-4.313 






46.73/20 


6.4 X 10-4 


3 - dep z 


1.0899 


0.467 


1.69 


-4.319 




0.028 


41.98/19 


0.020 


4 - dep L 


1.0652 


0.429 


1.70 


-4.386 


0.056 




30.19/19 


0.049 


5 - dep L,z 


1.0230 


0.369 


1.69 


-4.578 


0.109 


-0.066 


22.77/18 


0.200 






+0.040 


+0.18 


+0.110 


+0.019 


+0.024 






1(T errors 




-0.031 


-0.31 


-0.086 


-0.025 


-0.016 







expected ones. Figure 4 shows that the average distribu- 
tion function is a-symmetrical with a long tail at large 
Rjf values. 

We then divided the Lx-z-Rjc space into 23 indepen- 
dent bins, which were chosen in order to look for possible 
dependences of the P ( Rx \Lx, z) shapeonthe luminosity 
and redshift (see Figure 3 for a representation of the Lx- 
z regions) while ensuring that about, at least, 10 AGN 
were observed in each bin (see Figure 5). 

After several trials, we found that the probability 
distribution function distribution of Rx is, indeed, a- 
symmetrical, showing a maximum at Rx=Ro, where the 
median is located. At Rx larger than Rq the distribu- 
tion is fairly well represented by a Lorentz function hav- 
ing width 7r, which provides a shallow decline at large 
Rx values, while at Rx smaller than Rq the exponent 
2 of the Lorentz function is better substituted by an ex- 
ponent 4 which gives a wider shoulder at Rx^Ro and 
then a steep decline at even lower Rx- The width of 
the left {Rx<Ro) part of the distribution is controlled 
by the 7; parameter. Therefore the probability distribu- 
tion function P{Rx\Lx, z) is expressed by the following 
formula: 



P{Rx) = { 



N_ 



A N 



l+(^ 



x-«o(^x.^)' 



1 



{Rx < Ro) 



{Rx > Ro), 



(2) 



where, in order to obtain a continuous function at 
Rx=Ra results A = \J "ir/^i, and the parameter N 
is constrained by the probability normalization require- 
ment: / P{Rx\Lx,z)dRx = 1. For 0.3<.2< 3.0 and 
42.2<Lo5'Lx<47.0, we allowed to vary, as a function of 
Lx and z, the position of the maximum (median) i?o of 
the distribution, according to the following formula: 

i?o = Rc [aL{LogLx - 44) + 1] [a,{z - 0.5) + 1] . (3) 

At redshifts and luminosities outside these ranges, Rq 
was kept constant, equal to the values assumed at the 
limits of the ranges. 

In Table 2 the results of the fits carried out using 
this parameterization arc reported. Confidence regions 
of each parameter were obtained by computing at a 
number of values around the best-fit solution, while leav- 
ing the other parameters free to float (see Lampton et al. 
1976). The 68% confidence regions quoted correspond to 

The solution without dependences on both the lumi- 
nosity and the redshift is rejected by the test, while 



the solutions either depending only on luminosity or red- 
shift provide barely acceptable fits to the data. A fairly 
good fit to the data (20% x^ probability) is instead pro- 
vided by the solution ^5 where both a dependence on 
the luminosity and the redshift is allowed. However, it 
should be noted that the parameter of the redshift de- 
pendence, az, is different from zero only at 3cr confidence 
level. 

The data and the shape of the best fit probabil- 
ity distribution function are shown in Figure 5, while the 
corresponding dependences of i?o on Lx and z are shown 
in Figure 6. In Figure 7 we show the shape of the best 
fit #5 probability distribution function in different bins 
of X-ray luminosity and redshift, with evidence (the con- 
tinuous lines) to the part which is actually constrained 
by the data. 



0.1 
0.01 

0.1 

0.01 
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Rx = log Lr/Lx 

Fig. 5. — Probability distribution function of Rx as a function 
of Lx and 2. The dashed line shows our best fit solution #5 (see 
Table 2). The limits in the Lx-2 space of the seven regions used 
are shown in Figure 3. 
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Fig. 6. — Dependence of the position of the median Rq of the 
probability distribution function on hx and z as computed in our 
best fit solution #5 (red continuous line; see text and Table 2). 

3.3. Dependence on Nh 

We did not find any significant dependence of the 
P{Rx\Lx , z) distribution on the N// column densities. 
The sample was divided into absorbed (N// >10^^ cm~^) 
and un-absorbed AGN. For both sub-samples the test 
on the best fit solution #5 (even with a different sam- 
pling in order to always observe at least 10 objects in 
each bin) provided probabilities larger than 10% (up to 
70%), and was then not able to reject our best fit distri- 
bution. This result is in agreement with the analysis of 
the X-ray absorption properties of the faint radio sources 

z = 0.3 z = 0.6 z=1.0 z=1.5 z = 2.0 z = 2.5 




-4-2 -6-4-S -6-4-2 -6-4-2 -6-4-2 -6-4-2 
R, = log {l^/L^) 

Fig. 7. — Probability distribution function in different bins of hx 
and z, as computed in our best fit solution #5 (see Table 2). The 
continuous line shows the range of Rx where the fit is constrained 
by the data, while the dashed line shows where the distribution is 
extrapolated. 
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Fig. 8. — 1.4 GHz luminosity as a function of the intrinsic 2-10 
keV luminosity of all AGN contained in our sample (only the AGN 
with logLx>42 have been used in our analysis). Radio detection 
are shown by green dots, while radio upper limits are shown by 
arrows. The dependence, according to fit #5, of the position of 
the median Rq of the P{Rx\Lx , z) distribution (see eq. 3) at 
redshift and 3 is shown by red continuos lines. For comparison the 
relations derived by Brinkmann et al. (2000; orange long dashed), 
Panessa et al. (2007; blue dotted), and Ballantyne (2009; magenta 
dashed) are shown. The relations derived by the measure of the 
AGN fundamental plane by Merloni et al. (2003) and Li et al. 
(2008) are shown (green dot-dashed and cyan dot-dashed lines, 
respectively) assuming a fixed BH mass of 10* Mq (see text). 

of the CDFS by Tozzi et al. (2009). 

3.4. Sanity checks 
3.4.1. The Lr-Lx relation 

We compared our measure of the probability distribu- 
tion function of Rx with previous measures on the re- 
lationship between the AGN radio and X-ray luminosi- 
ties. In Figure 8 we show the 1.4 GHz luminosity as a 
function of the intrinsic 2-10 keV luminosity for all the 
sources of our sample. Our analysis differs from many 
previous studies where single general log-log linear re- 
lations were derived. In these studies, even when the 
presence of censored data, or data with errors in both 
axes (see e.g. La Franca et al. 1995), were taken into 
account, it was assumed the presence of a symmetrical 
(usually Gaussian) distribution of the deviations from 
the best fit relation, which were attributed to an intrin- 
sic scatter. Our method, instead, by taking into account 
all the censored data, allows to measure the shape of the 
distribution of the intrinsic scatter and its possible de- 
pendences on other variables (such as the luminosity and 
redshift in our case) . In Figure 8 we compare our best fit 
solution #5 with the relations derived by other authors 
such as Brinkmann et al. (2000) in the soft 0.5-2 keV 
band^, Merloni et al. (2003), Panessa et al. (2007), Li 
et al. (2008) and Ballantyne (2009; for the radio quiet 
AGN) in the 2-10 keV band. The relations by Merloni et 

^ An X-ray spectral index a=0.7 was assumed to convert the 
0.5-2 keV luminosities into 2-10 keV luminosities. 
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al. (2003) and Li et aL (2008) were derived from their 
measure of the AGN fundamental plane (i.e. including 
also a dependence on the BH mass), assuming a fixed 
BH mass of 10^ Mq. 

As our measured distribution function is a-symmctrical 
and depends on the redshift, the comparison is not 
straightforward (in Figure 8 we plot the position of the 
median Rq of the distribution at redshifts and 3). It 
results that our best fit solution #5 of the dependence 
of the 1.4 GHz luminosity on the 2-10 keV luminosity is 
flatter than obtained in previous works. Our fitted re- 
lation corresponds to a power-law L/jocL^ with indexes 
a=0.48 and a=0.58 at redshifts and 3, respectively. 
The relations fitted by the other authors have instead a 
wide range of power-law indexes (0.6< a < 1-5), system- 
atically steeper than our result. This difference is partly 
caused by our new method used to measure the L^j-Lx 
relation but it is also caused by our introduction (and 
measure) of a dependence on the redshift of the average 
Hx values. In fact, as shown in Figure 8, our solution 
T^4, obtained without the inclusion of a redshift depen- 
dence, results in a steeper slope, having a power-law in- 
dex a = 0.75. Figure 8 also shows that the fundamental 
plane by Merloni et al. (2003) has a power-law index 
(a=0.60) close to our best fit estimate (a=0.48-0.58), 
and for BH masses of 10^ Mq is located between our 
best fit relations at redshift and 3. In the fundamental 
plane estimate by Merloni et al. (2003) Kx has also a 
dependence on the BH mass of the type oc0.781og(MBH)- 
This implies that our estimates at redshift and 3 are 
roughly similar to the fundamental plane measures for 
BH masses of 10^ and 10^ M0, respectively. Therefore, 
as in our (and all fiux limited) samples, high redshift 
AGN are on average more luminous, and thus probably 
host on average more massive BH, we can infer that our 
measure of the increase with the redshift of the median 
of the Kx distribution is qualitatively in agreement with 
the AGN fundamental plan measures where an increase 
of the average Kx with the BH masses is observed. 

3.4.2. The I.4 GHz luminosity function and counts 

We then verified if the measured probability distribu- 
tion function of Kx, once convolved with the 2-10 keV 
LF, properly reproduces (taking into account the uncer- 
tainties) independent previous measures of the 1.4 GHz 
LF, $fl;(Li.4,z), and integral counts, A^(> S). We used 
the 2-10 keV LF, ^x{Lx,z), as measured by La Franca 
et al. (2005), and modified by allowing a steep exponen- 
tial decline of the AGN density at redshifts larger than 
z = 2.7, as measured by Brusa et al. (2009a). As dis- 
cussed in La Franca et al. (2005), the density of Comp- 
ton thick AGN with 24<logN//<26 cm~^ was assumed 
to be equal to the density of the Compton thin AGN with 
22<logN//<24 cm^^ (this assumption resulted to prop- 
erly reproduce the cosmic X-ray background). The X-ray 
LF can be converted into a radio LF by the formula: 

$i?(ii.4,z) = J P{Rx\Lx,z)<^xiLx,z)d\ogLx. (4) 

The X-ray LF was integrated starting from an X-ray lu- 
minosity logLx=41 erg/s. The predicted radio LF is 
shown in Figure 9 and compared with the 1.4 GHz ra- 
dio LF measured by Smolcic et al. (2009). As discussed 
in Smolcic et al. (2009), this radio LF includes (like in 
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Fig. 9. — Predicted 1.4 GHz radio LF according to our fit #5, at 
redshifts 0.2, 0.4, 0.7 and 1.0 (red continuous lines). The dashed 
red lines show the part of the LF reproduced by the P(Rjf ) distri- 
bution with Rx>-4. The black long dashed lines show the radio 
luminosity function as estimated by Smolcic et al. (2009). 

our measure of the distribution of Kx) mostly the FRI 
sources. Although also affected by the uncertainties on 
the measure of the X-ray LF by La Franca et al. (2005), 
our best fit #5 of the probability distribution function of 
Kx provides a fairly good reproduction of the FRI radio 
luminosity function. 

As a result of this computation, our estimate of the 
distribution of Kx allows to predict the AGN FRI ra- 
dio LF at luminosities smaller than Li.4=10'^^ W Hz~^, 
never probed before. However at these low radio lumi- 
nosities the AGN could be overwhelmed by the radio 
emission of the hosting galaxy because during periods of 
strong star formation activity the supernovae remnants 
accelerate cosmic rays which radiate synchrotron emis- 
sion in local magnetic fields. According to Ranalli et 
al. (2003), strong star forming galaxies (~10^ Mo/yr) 
have about 10^^ erg/s 2-10 keV luminosities, and, in gen- 
eral, in all star forming galaxies a linear relation between 
the 1.4 GHz radio and the 2-10 keV luminosites is ob- 
served, which corresponds to a value Kx— —2.0. At 
face value this would imply that for the lowest lumi- 
nosity AGN (logLjf— 42-43 erg/s) most of our measure 
of the P(Rx) distribution (which spans in the range 
— 7<Rx< —1) could be contaminated if the hosting 
galaxies are undergoing a strong star formation activ- 
ity. However, as shown in Figure 7, at low Kx values the 
P(Rx) distribution is mostly measured from AGN with 
43<logLx<45 erg/s, and redshift lower than z~0.5. At 
these low redshifts it is very unlikely for galaxies to har- 
bor a star formation stronger than ^10 Mq/yi (Elbaz et 
al. 2007; Noeske et al. 2007), which roughly corresponds 
to an emission in the 2-10 keV band of logLx~41 erg/s 
(Ranalli et al. 2003). Therefore, for example, in a low 
redshift AGN with logLx~44 erg/s, having an hosting 
galaxy whose star formation emits a 2-10 keV luminos- 
ity of logLx^41 erg/s, the corresponding radio emission 
(having an intrinsic Rx— -2.0) can contaminate our mea- 
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Fig. 10. — Euclidean differential radio counts at 1.4 GHz. The 
black continuous line shows the expected counts according to the 
convolution of the P(Rx) distribution of fit # 5 (see Table 2) with 
the X-ray LF of La Franca et al. (2005). The green dashed line 
shows the radio counts produced by the FRII sources according 
to Wilman et al. (2008). The red continuous line is the sum of 
our predicted counts (the black line) with the counts of the FRII 
sources (green dashed line). Radio counts measures are shown 
by black open circles (FIRST; White et al. 1997), magenta filled 
squares (Bondi et al. 2008) and orange filled circles (Owen et al. 
2008). Blue open squares and cyan open triangles show the AGN 
counts estimated by Seymour et al. (2008) and Padovani et al. 
(2009), respectively. 



sure of the Kx distribution only for Kx values smaller 
than -5. We can then conclude that only at the lowest 
Kx values (< —5) our measure of the P(Rx) distribu- 
tion is potentially affected by a contamination from the 
radio emission due to the star formation activity of the 
AGN hosting galaxies. 

We used the above computed radio LF (eq. 4) to derive 
the expected integral counts N(>S) from the following 
equation: 



$(^1.4, 2;)cilogL 



1.4, 



(5) 

where k{z) = (1 + z)"""^ (with oc and or = 

0.5) is the radio k-correction, di{z) the luminosity dis- 
tance, and the integral counts are measured in sr^^ units. 
The predicted radio counts are shown in Figure 10. As 
our measured distribution of Rx represents only the FRI 
population (see discussion in sect. 2), in order to re- 
produce the total radio counts, the contribution derived 
from the LF of the FRII population (as measured by 
Wilman et al. 2008) was added. The reproduced counts 
are in good agreement with the observations. At 1.4 GHz 
fluxes below 1 mJy the euclidean radio counts flatten, due 
to the appearance of the population of the star forming 
galaxies. In this context, it is matter of discussion which 
is the fraction of the AGN at these fluxes. Our results 
agree with recent estimates of the AGN contributions 
to the sub-milliJansky radio counts from Seymour et al. 

(2008) and Padovani et al. (2009). Similar results (at 
these fluxes) have also been obtained from semi-empirical 
simulations of the extragalactic radio counts by Jarvis & 
Rawlings (2004), Wilman et al. (2008) and Ballantyne 

(2009) . 
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Fig. 11. — Kinetic luminosity as function of the radio luminosity 
as estimated by Willott et al. (1999; magenta lines), Best et al. 
(2006; blue Unes), Merloni & Heinz (2007; red dotted Une) and 
Birzan et al. (2008; green dotted line). In our work we used the 
combination of the relations from Best et al. (2006) and Willott 
et al. (1999) at low and high luminosities, respectively (the two 
continuous lines). 

4. THE KINETIC LUMINOSITY FUNCTION 

Once the radio LF of the FRI sources is measured 
(Equation 4 and Figure 9) , in order to derive the kinetic 
mechanical LF and its evolution, we should convolve the 
1.4 GHz LF with a relation which converts the radio lu- 
minosity Li.4 into a mechanical power \jk- In the last 
decade, several authors have worked on the estimate of 
this relation (Willott et al. 1999; Birzan et al. 2004, 
2008; Best et al. 2006; Heinz et al. 2007; Merloni & 
Heinz 2008). Following the discussion of Cattaneo and 
Best (2009) we used two different estimates which are 
representative of two different luminosity regimes. At 
high radio luminosities (Li.4>10^^ W Hz~^), Willott et 
al. (1999) used the minimum entropy density that the 
plasma radio lobes must have in order to emit the ob- 
served synchrotron radiation and obtained (see Figure 
11): 

0.85 



Lk = 1.4 X 10' 



37 



( Li. 



4 GHz 



VlO^sWHz" 



W. 



(6) 



A second approach is to infer Lk from the mechani- 
cal work that the lobes do on the surrounding hot gas. 
The expanding lobes of relativistic synchrotron-emitting 
plasma open cavities (of volume V) in the ambient ther- 
mal X-ray emitting plasma. The minimum work in in- 
flating these cavities is done for reversible (quasi-static) 
inflation and equals pV ^ where p is the pressure of the 
ambient gas. Best et al. (2006) derived a relation be- 
tween radio and mechanical luminosity based upon this 
estimate for the energy associated with these cavities, 
combined with an estimate of the cavity ages from the 
buoyancy timescale (from Birzan et al. 2004). Com- 
paring the mechanical luminosities of 19 nearby radio 
sources that have associated X-ray cavities with their 
1.4 GHz monochromatic radio luminosities leads to a re- 
lation 



Lk = 1.2 X 10' 



,37 



4 GHz 



0.40 



W, 



(7) 



Vl025WHz' 

W Hz^^) and is close to the estimate by Birzan et al. 



which is better suited for low luminosities (L1.4 
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Fig. 12. — Bolometric radiative (green continuous line) and Iti- 
netic AGN LF. The magenta and blue lines show the kinetic LF 
predicted using the relations linking the radio and kinetic lumi- 
nosities of Willott et al. (1999) and Best et al. (2006) respectively. 
The red line is the result obtained using the combination of the 
relations of Willott et al. (1999) and Best et al. (2006) at high 
and low luminosities respectively (see text and Figure 11). 

(2008; see Figure 11). Using a similar method, Merloni 
& Heinz (2008) obtained at high radio luminosities, a 
relation similar to the one by Willott et al. (1999) but 
with about 0.5-1 dex higher kinetic luminosities (see Fig- 
ure 11). 

In this work (as in Cattaneo and Best 2009) wc used 
Equation 6 at high luminosities (Li.4>10^^ W Hz^^) and 
Equation 7 at lower luminosities. Using these relations, 
once derived the radio LF from the X-ray LF and the Kx 
distribution according to Equation 4, we can estimate the 
kinetic LF, ^k{Lk, z), by the formula 



dNjLK.z) 
dVdlogLK 



$1.4(il.4(iK),z) 



dlogLi 
dlogLx 



(8) 

while the bolometric radiative LF, ^radiLs, z), is derived 
from the X-ray LF, via the bolometric correction Lg = 
Kx{,Lx)Lx, by 



^rad{LB,z) 



dN{LB,z) 
dVdlogLB 



^x{Lx{Lb):z) 



dWogLx 



dlogLB ' 
(9) 

where for Kx (L^ ) we used the relation from Marconi et 
al. (2004). Both the, so derived, kinetic and radiative 
(bolometric) LF are shown in Figure 12, while in Figure 
13 the kinetic power density, Lik^k(J^k ,z), as function 
of the kinetic luminosity is shown. 
By integrating in luminosity the above derived kinetic 
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Fig. 13. — Kinetic power density, hK^ni^K as function of 
the kinetic luminosity according to the fit solution #5 (red contin- 
uos line). The results by Merloni & Heinz (2008) are shown by a 
cyan shaded area. 

and radiative (bolometric) LF it is possible to estimate 
the dependence of the AGN mechanical and radiative 
power per unit cosmic volume as a function of redshift. 



^k{z) = J LKiLiA)^iA{LiA,z)dlogLiA, (10) 
^rad{z)^ j Lboi{Lx)^x{Lx,z)d\ogLx, (11) 



and 



which are shown in Figure 14. As both the above equa- 
tions depends on the X-ray luminosity function (see eq. 
4) the integral lower limit was logLx=41 erg/s which 
corresponds to a lower limit of of logLB~35.5 W for eq. 
11 and, according to the Rx distribution and the rela- 
tion between L1.4 and L^ (eq. 7) to a lower limit of 
logL;f~34.5 W for eq. 10. The resulting kinetic power 
density would be similar if, at all luminosities, the con- 
version of the radio into kinetic luminosity by Best et al. 
(2006), eq. 7, would be used, while a factor 4-6 lower val- 
ues would be obtained using the relation from Willott et 
al. (1999), eq. 6, only. As shown in Figure 12 this differ- 
ence is caused by the steep drop off of the 1.4 GHz radio 
LF at luminosities higher than LogLfl^25 W (see Fig- 
ure 9); as a consequence the resulting power density (the 
integrated kinetic LF) depends mainly on which conver- 
sion of the radio luminosity into kinetic power is used at 
low luminosities. 



Radio loudness and the AGN feedback 



11 



this work 




: ^km(^) cav + synchro 

: "k,„(z) (Rx>-'') 



I I I I I I I I I I I 1 I 

1 2 3 4 5 6 

redshift 

Fig. 14. — Bottom. Power density as a function of redshift. 
The kinetic power density derived from the best fit solution #5 
is shown by a continuous red line. The red dashed line shows 
the result if the more radio loud objects (Rjsf >-4) only are used 
(see text and Figure 9). The magenta and orange dashed lines 
show the kinetic luminosity density as estimated by Kording et 
al. (2008) and Cattaneo & Best (2009), respectively. The cyan 
and blue areas show the estimates by Merloni & Heinz (2008) and 
Smolcic et al. (2010), respectively. The radiative and kinetic power 
densities, as used in the model of Croton et al. (2006), are shown 
by continuous and dashed black lines, respectively. Top. Ratio 
between the kinetic and radiative power densities. The red and 
black continuous lines are the results by our fit #5 and Croton et 
al. (2006), respectively. 

5. DISCUSSION 

The probabihty distribution function of Rx, 
P{Rx\Lx,z) estimated in Section 3 depends on 
both luminosity and redshift: the average Rx increases 
with decreasing luminosity and (possibly) increasing 
redshift (see best fit #5 in Table 2 and Figure 6). 
The observed increase of the average Rx value with 
decreasing luminosity is similar to previous results and 
models where (in analogy with X-ray binaries) low 
luminosity AGN are expected to be more likely radio 
loud (see e.g. Merloni & Heinz 2008; Kording, Jester 
& Fender 2008). 

The knowledge of the Rx distribution, once convolved 
with the X-ray LF and the relations between the kinetic 
and radio luminosity, allows to estimate the kinetic LF 
and its evolution. At luminosities higher than the break 
of the bolometric LF (L^^IO^^ W) the kinetic LF resuhs 
to be more than two orders of magnitude smaller than the 
bolometric LF (see Figure 12), while at lower luminosities 
the relevance of the kinetic LF increases, reaching values 
comparable to the bolometric one at L/f-lO^^ W, which 

* As discussed in sect. 2 and 3, these trends do not include 
the most luminous (radio loud; mostly FRIl) sources. As far as 
the more radio loud population is concerned, it has been observed 
that the fraction of more radio loud AGN increases with increasing 
optical (or X-ray) luminosity and decreasing redshift (Miller et al. 
1990; Visnovsky et al. 1992; Padovani 1993; La Franca et al. 1994, 
Goldschmidt et al. 1999, Ivezic et al. 2002; Cirasuolo et al. 2003; 
Jiang et al. 2007). 



roughly corresponds to the minimum luminosity experi- 
mentally probed by the X-ray LF (logLx'^42 erg/s; see 
e.g. La Franca et al. 2005). The kinetic LF shows a 
maximum in the range 10'^^<Lk<W'^^ W, where most 
('--^90%) of the kinetic power density (shown in Figure 
13) is produced. 

In Figure 14 we show the kinetic and radiative power 
density as a function of the redshift. We also show the 
kinetic power density corresponding to the more radio 
emitting AGN having Rx larger than -4, and then corre- 
sponding to the population which is typically represented 
by the radio LF (see in Figure 9 the comparison of the 
reproduced radio LF with Rx>-4 with the radio LF of 
FRI sources as measured by Smolcic et al. 2009). It is 
then possible to see that the kinetic power density could 
be underestimated by up to a factor of about two, if the 
radio LF alone is used, without taking into account the 
low radio luminosity AGN population. 

Our estimates are in qualitative agreement with the 
trends of the radiative and kinetic power density with 
redshift used by Croton et al (2006) at z > 0.5 (see Fig- 
ure 14). However, at lower redshifts we find a sharp 
(a factor of five) decrease of both radiative and kinetic 
power densities from z ~0.5 (i.e. about 5 billion years 
ago) to z = 0. This result is quite robust, and comes 
from the strong negative evolution of the AGN LF from 
z 2 down to 2 = 0, which has been observed since 
the first studies of the QSO evolution in the optical (see 
e.g. Marshall et al. 1983; Groom et al. 2009 for re- 
cent results), and measured by many other authors in 
the hard X-rays (e.g. Ueda et al. 2003; La Franca et al. 
2005; Hasinger 2008). Conversely, Groton et al. (2006) 
assume an almost continuous increase of both the kinetic 
and radiative powers, due to the assumption that both 
phenomena are related to an almost constant accretion 
onto the SMBH. Under this assumption Groton et al. 
(2006) overestimate both the AGN radiative and kinetic 
power densities at low (z<0.5) redshift, allowing only for 
a much shallower decrease of the kinetic feedback (a fac- 
tor of 30%) and of the AGN radiative power. 

In Figure 14 (top) we show the ratio of the ilx/^rad 
plotted as function of redshift. According to our best fit 
#5, at z>0.5 the kinetic power density is ^ 5% of the 
radiative density, and increases up to about 30% at de- 
creasing redshifts. This increase at low redshifts of the 
ratio of the kinetic power to the radiative power density 
could help in modeling the quenching of the star forma- 
tion at low redshift. In Groton et al. (2006) a milder 
increase is assumed, which should be attributed to the, 
above discussed, overestimate of the low redshift AGN 
radiative density. 

Many previous results on the AGN kinetic LF are 
based on the convolution of some relations between the 
kinetic and radio power with direct measures of the AGN 
radio LF (e.g. Shankar et al. 2008; Merloni & Heinz 
2008; Kording, Jester & Fender 2008; Cattaneo & Best 
2009; Smolcic et al. 2009). As already discussed in the 
introduction, the measure of the Rx distribution is use- 
ful in order to allow a detailed implementation of the 
AGN feedback within the galaxy formation and evolu- 
tion models because it gives the opportunity to predict 
the radio luminosity (and thus feedback) of each AGN as 
a function of its luminosity (accretion rate) and redshift. 

According to our best fit #5, in the redshift range 
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0.5<z<3 the integrated kinetic power density is 
2x 10"^^ W Mpc^^ (see Figure 14). This is in rough agree- 
ment with the previous estimates by Merloni & Heinz 
(2008), Kording, Jester and Fender (2008) and Smolcic 
et al. (2009). At lower redshift {z<0.5) we observe a drop 
by a factor of five, similar to what observed by Kording, 
Jester and Fender (2008), and in agreement, within the 
uncertainties, with Merloni & Heinz (2008). On the con- 
trary our results are, at any redshift, 2-8 times greater 
than that reported by Cattaneo & Best (2009). 

Merloni & Heinz (2008) found that their kinetic LF 
roughly corresponds to a constant overall efficiency in 
converting the accreted mass energy into kinetic power 
ekin—3 — 5 X 10~^ (where Lk = €kinmc2). Their results 
are similar to ours where, according to equations 10 and 
11, on average, we have €kin—{^K/^R)^R — 5 x 10"^, 
as we measure VLk/^b. — 0.05 (see the VIk/^r ratio as 
a function of redshift in Figure 14, top), and assuming a 
radiative efficiency Ci? = 0.1 (Marconi et al. 2004). How- 
ever, we observe an increase of the / ratio (i.e. of 
the kinetic efficiency) up to a value 0.3 at decreasing 
redshifts, which Merloni & Heinz (2008) observe in the 
most massive objects only (> 10* — lO^M©; see e.g. their 
Figure 13). It should also be noted that, although the in- 
tegrated kinetic power density of Merloni & Heinz (2008) 
is in agreement with our estimate, their kinetic LF is sim- 
ilar to our measure only at L^^IO'^^ W (see Figure 13), 
while it is definitely larger (by about an order of mag- 
nitude) at higher kinetic hnninosities. Once integrated 
in luminosity, the computed power densities are similar 
(at 2;<3) because our low luminosity limit (Lis:=10^^ W) 
is significantly lower than that used by Merloni & Heinz 
(2008; hK=lQ^^ W). 

Shankar et al. (2008), found that the ratio of the ki- 
netic to bolometric luminosity, defined as Qk = Lk/Lb, 
is constant and equal to gk = 0.10 with a scatter of 
a = 0.38. According to equations 10 and 11, gk cor- 
responds roughly to the ratio D,K/^rad (plotted as func- 
tion of redshift in Figure 14, top), which, as discussed 
above, levels at ~0.05 at z > 1 while increases up to 0.3 
at z = 0. 

6. CONCLUSIONS 

We used a sample of more than 1600 X-ray selected 
AGN observed at 1.4 GHz to measure the probability 
distribution function, P{Rx\Lx, z), of the ratio Rx of 
the radio to intrinsic X-ray luminosity, as a function of 
the X-ray luminosity and redshift. 

The knowledge of the P(Rx\Lx, z) distribution is nec- 
essary to estimate the AGN kinetic (radio) feedback into 
the hosting galaxies by allowing to couple it with the 
luminous, accreting, phases of the AGN activity. 

The average value of Rx increases with decreasing X- 
ray luminosities and (possibly) increasing redshift. At 
variance, we did not find a statistical significant differ- 
ence between the radio properties of the X-ray absorbed 
(Nh>1022 cm-2) and un-absorbed AGN. 

We were able to better measure the densities of the 
more radio quiet (Rx<-4) AGN which resulted to be 
responsible of about half of the derived kinetic power 
density. 

According to our analysis the value of the kinetic en- 
ergy density is in qualitative agreement with the last 
generation galaxy evolution scenarios, where radio mode 



AGN feedback is invoked to quench the star formation 
in galaxies and slow down the cooling flows in galaxy 
clusters. However at redshifts below 0.5, similarly to 
what observed by Kording, Jester and Fender (2008), we 
find a sharp (about a factor of five) decrease of the ki- 
netic energy density, which is strictly related the AGN 
density evolution, but which is not included in many of 
the galaxy/ AGN formation and evolution models where, 
instead, the radio mode feedback is assumed to contin- 
uously increase (or only smoothly decrease) at low red- 
shift. 
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machine readable format. We thanks the anonymous ref- 
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